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Abstract 


Simulation results obtained by using FUN2D for robust airfoil shape optimization in tran- 
sonic* viscous flow are included to show the potential of the profile optimization method for 
generating fairly smooth optimal airfoils with no off-design performance degradation. 
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chord length of airfoil 
drag coefficient. 

gradient of c d with respect to D 
derivative of c<i with respect to a 
lift coefficient 

gradient of c\ with respect to D 

derivative of c/ with respect to a 

target lift coefficient 

design vector 

mean of random variable 

feasible set for the design vector D 

free-stream Mach number 

number of design variables 

probability density function of Mach number 

number of design conditions 

coordinates of points on plane 

angle of attack 

minimum rate of drag reduction at all design conditions 

change in design vector 

scalars defining the trust region 

variance of random variable 

a given Mach range 

inner product in Euclidean space 


Subscripts and Superscripts 

i index for design condition 

j index for component of design vector 

k index for iteration or iterate 


1 Introduction 

This paper is an extension of reference 1, where the profile optimization method was pro- 
posed to resolve off-design performance degradation problems encountered by multipoint 
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optimization methods (ref. 2). Numerical simulation results given in reference 1, though 
very encouraging, were not representative of realistic airfoil shape optimization problems. 
To strengthen the simulation results in reference 1, we apply the profile optimization method 
to two cases of airfoil shape optimization in transonic viscous flow over a range of Mach 
numbers. 

The first case is the redesign of the RAE2822 airfoil (ref. 2) over the Mach range from 0.68 
to 0.76, with the target lift at 0.733 and Reynolds number 2.7 x 10 6 . The second case is the 
redesign of Whitcomb’s integral supercritical airfoil (ref. 3) over the Mach range from 0.68 
to 0.77, with the target lift at 0.7 and Reynolds number 2.7 x 10 6 . Viscous flow analysis and 
the corresponding gradient calculation are based on FUN2D (ref. 4), which can perform fully 
turbulent Navier-Stokes flow analysis and the corresponding discrete adjoint analysis with 
unstructured grids. Two spar location thickness constraints and the maximum thickness 
constraint are enforced for both airfoil shape 1 optimization cases. Airfoils are parameterized 
by 35 B-spline control points. Design variables are angles of attack at specified design 
conditions and the ^-coordinates of 35 B-spline control points. 

In both cases, the simulation results show that the profile optimization method can 
generate fairly realistic optimal airfoil shapes in transonic viscous flow without severe off- 
design performance degradation, when 4 design conditions and 35 geometric* design variables 
are used. 

Because the purpose of this paper is to provide additional simulation results for the 
profile^ optimization method given in reference 1, readers should consult reference 1 for 
relevant references and technical details. 

The paper is organized as follows. In section 2, we review the profile optimization 
method. Section 3 includes the simulation results for airfoil shape optimization in transonic 
viscous flow. Concluding remarks are given in the final section. 

2 Profile Optimization Method 

In this section, we briefly review the profile optimization method proposed in reference 
1. The profile optimization method is intended to solve the following robust optimization 
problem: 

A ( E(c ■ t> •'' <e ^‘ , ) (I) 


DZT and c(D,a(A/),Jlf)=cf(M) for M € SI (2) 


Here cf(M) is the target lift requirement for Mach number M, T is a given feasible set for 
geometric design variables (that could be defined by geometry constraints such as thickness 
constraints), and a (A/) is the angle of attack corresponding to Af. The drag and lift 
coefficients are Cd and C/, respectively. The mean and variance of Cd are defined as 


E(c d ) 


-j 

Ju 


c d (D,a(M),M) ■ p(M) dM, 


<r 2 (c d ) = j ^ [r.d(D,a(M), M) - E(c d )} 2 p(M) dM , J 


( 3 ) 


where p(M) is a probability density function of M and 0 is a given Mach range (such as 
from M = 0.68 to M = 0.77). 

The robust optimization model in equation (1) addresses some important issues in aero- 
dynamic shape optimization. For example, to avoid off-design performance degradation, one 
can reduce a(cd) as much as possible. Note that if rr(cd) = 0, the corresponding solution will 
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have the same (perhaps poor) performance over the given Mach range. Therefore, equation 
(1) provides a tool for trade-offs between average performance improvement and perfor- 
mance fluctuations over the Mach range. However, due to the high computational cost for 
solving Navier-Stokes equations, it is impossible to have reasonable estimates of E(cd) and 
rr(crf). Consequently, it is not practical to solve equation (1) using current computational 
fluid dynamics tools. 

One alternative way to avoid off-design performance degradation is to find a descent 
direction that could reduce the drag simultaneously and proportionally over the given Mach 
range while keeping the lift at the target value. The profile optimization method was de- 
signed to find such a descent direction with limited information on lift and drag over the 
given Mach range. The innovative feature of the profile optimization method is to adaptively 
change the objective function from iteration to iteration to achieve simultaneous and pro- 
portional drag reduction over the given Mach range. In contrast to methods that minimize 
one aggregate objective function to find a Pareto optimal solution to a multiobjective opti- 
mization problem, the profile optimization method does not solve any optimization problem 
with one objective function; instead, it looks for a Pareto optimal solution that has the least 
chance for severe off-design performance degradation. 

Profile Optimization Method. Let D° be a given initial design vector, let M \, . . . , M r 
be a set of design points over the given Mach range, and k — 0. Construct a sequence of 
design vectors as follows: 

1. Compute feasible angles of attack. Find . . . , a r .* such that 

ci{D k ,ai, k ,Mi) = c* t i for 1 < i < r. 

2. Solve a trust region subproblem. Let c-d,i,k and be the linear approximations of 
the drag and lift at (D*,**,-,*): 

«W(A0, AaO = c,(D k ,a,, k , M«) + (^, AD^ + ^Aa ; , 

c d j, k (AD,A ai ) = Cd(D k ,a it k,Mi) + AD^ + ^Aa i; 

where the derivatives are evaluated at (D k ,{**.*, M*). Consider the following trust 
region subproblem: 

min r / subject to D k + AD € T , (4) 

A/XAa, 

&i,kpk 5^ AD i ^ &i,kPk for 1 ^ i ^ w, 

-a* < Aa 7 < a*. for 1 < i < r, 

ci,i,k(AD, A a<) = cl for 1 < i < r, 

Cd,i.k(AD, Aa,) < (1 - r)) ■c d (D k ,a i , k ,M i ) for 1 < i < r, 

where (5, ^ > 0 and a*. > 0 are scalars that determine the trust region, and D k 4- AD € 
T means that the airfoil corresponding to ( D k + AD) satisfies all the geometric con- 
straints. (For our simulation runs, Si ^ is approximately the thickness of the airfoil 
at the x-coordinate of the ith control point, a* = 1, and D k -f AD e T means that 
the airfoil satisfies two thickness constraints at the spar locations and the maximum 
thickness constraint.) Determine the smallest pk > 0 such that equation (4) is feasi- 
ble and the least norm solution (AD k , Aai,*, . . . , Aa r ^) of equation (4) satisfies the 
following condition: 

Cd,i,k{AD k , Aati'k) < (1 - lm\n)c d {D k ,a i<k ,Mi) for 1 < i < r. 
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3. Generate the now iterate. Let <*>.*+ 1 = a;.*. + An,-,*, for 1 < i < r and D /, + 1 — 
Z? A + AZ?*. 

4. Start a new iteration. Update A- by k - hi and go bark to step 1. 


3 Numerical Simulation Results 

We apply the profile optimization method in two cases of airfoil shape design optimization 
in transonic viscous flow. The first case is the redesign of the RAE2822 airfoil over the Mach 
range from 0.G8 to 0.76. We adopt the four design conditions used in Drela’s simulation 
(ref. 2): M — 0.68,0-71,0.74, and 0.76, with the target lift at 0.733 and Reynolds number 
2.7 x 10°. The second case is the redesign of Whitcomb’s integral supercritical airfoil (ref. 
3) over the Mach range from 0.68 to 0.77. In this case, we use the following four design 
conditions: M - 0.68,0.71,0.74, and 0.77, with the target lift at 0.7 and Reynolds number 
2.7 x 10 (i . 

For a given airfoil, the lift and drag coefficients and their gradients an' calculated by 
solving fully turbulent Navier-Stokes equations and the corresponding discrete adjoint equa- 
tions using FUN2D (ref. 4). In both cases, the unstructured grid has 300 grid points on 
the airfoil and 32 grid points on the far field (which is placed at 20 chord lengths). The 
unstructured grid has a total of 18654 grid points, 55631 elements, and 37308 faces (see fig. 
1 for the grid near the RAE2822 airfoil). 

With the given grid, the flow solver and the adjoint solver are terminated when the 
2-norms of the residual of the density equations and the corresponding discrete adjoint 
equations are reduced by at least five orders of magnitude. 

Airfoils are parameterized by 35 B-spline control points, which is large enough to ac- 
commodate free- form geometry changes for both cases. Figure 2 shows parameterization of 
the RAE2822 airfoil and Whitcomb’s integral supercritical airfoil. 

The ^-coordinates of all the control points are fixed during optimization. Changes of 
the ^-coordinates of the five control points near the trailing edge are constrained to be 
the same so that the shape of the airfoil near the trailing edge is not too oscillatory. The 
same geometry constraint is applied to the //-coordinates of the three control points near 
the leading edge. The reason for such geometric constraints is that the optimizer may not 
be able to make a reasonable modification of the shape of the leading or trailing edge at 
very fine scales. All the //-coordinates of the 35 control points are used as shape design 
variables. In both cases, we impose a thickness constraint at spar locations x/c = 0.15 and 
x/c = 0.6 and at the maximum thickness location. For each iterate D k , we find angles of 
attack a I.*,.. . ,a r ,ft such that |r/(D*, a 7 ,*, AJ<) - < 0.001. 


3.1 Drag Reduction Over the Given Mach Range 

Figures 3 and 4 show the history of drag changes at the design conditions for cases 1 and 2, 
respectively. 

Figure 5 shows the postoptimization analysis of drag curve over the given Mach range 
for cases 1 and 2. The drag rise curves are constructed by using 10 equally spaced Mach 
numbers from M - 0.68 to M = 0.77, which include 4 design points and 6 off-design 
points. In both cases, the profile optimization method reduces the drag of the baseline over 
the given Mach range. Moreover, there is no off-design performance degradation of the 
optimized airfoils. 
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3.2 Airfoil Shapes 

The optimal airfoils generated by the profile optimization method are quite smooth and 
do not have shock bumps corresponding to design conditions (i.e., Mach numbers). This 
result is different from the numerical simulation results obtained earlier by Drela (ref. 2) on 
multipoint airfoil optimization when sinusoidal basis functions were used to parameterize 
the airfoil shape. Figure 6 shows the optimal airfoils generated by the profile optimization 
method for cases 1 and 2. 

3.3 Pressure Distributions 

Pressure distributions for the baseline and the optimal airfoil in both cases are plotted in 
figures 7 and 8. 


4 Concluding Remarks 

This paper documents two cases of using the profile optimization method to redesign airfoils 
over a range of Mach numbers in transonic viscous flow with airfoils parameterized by 35 
B-spline control points. FUN2D (ref. 4) is used to compute the lift and drag values and 
their gradients with respect to changes in airfoil shape and changes in angle of attack. All 
the ^/-coordinates of 35 B-spline control points are used by the profile optimization method 
to search for the true optimal solution. With thickness constraints at two spar locations 
and the maximum thickness constraint, the profile optimization method minimizes the drag 
at four design conditions while keeping the lift at the target value. The optimized airfoils 
have fairly realistic and smooth shapes. Postoptimization analysis shows that the optimized 
airfoils have no off-design performance degradation. 

The simulation results demonstrate that the profile optimization method has the poten- 
tial to be used for real-world airfoil shape optimization over a range of flight conditions. 
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Figure 1: Grid used to solve Navier-Stokes equations by FUN2D. 




Figure 2: RAE2822 airfoil and Whitcomb’s integral supercritical airfoil parameterized by 
35 cubic B-spline control points. 


















Figure 7: Pressure distributions on RAE2822 and the corresponding optimal airfoil at four 
design conditions with target lift at 0.733. 






Figure 8: Pressure distributions on Whitcomb’s integral supercritical airfoil and the corre- 
sponding optimal airfoil at four design conditions with target lift at 0.7. 
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